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Abstract 

Superheavy particles of masses ~ 10 13 — 10 14 GeV with lifetimes ~ 10 10 — 
10 22 years are very interesting, since their decays may account for the ultra-high 
energy (UHE) cosmic rays discovered beyond the Greisen-Zatsepin-Kuzmin cut- 
off energy E ~ 5 x 10 10 GeV. We show that the presence of such long-lived 
superheavy particles is a generic prediction of QCD-like S\J(N C ) gauge theories 
with Nf flavors of quarks and antiquarks and the large number of colors N c . We 
construct explicit models based on supersymmetric SU(iV c ) gauge theories and 
show that if the dynamical scale A ~ 10 13 — 10 14 GeV and N c = 6 — 10 the 
lightest composite baryons have the desired masses and lifetimes to explain the 
UHE cosmic rays. Interesting is that in these models the gaugino condensation 
necessarily occurs and hence these models may play a role of so-called hidden 
sector for supersymmetry breaking in supergravity. 



1 Introduction 



In Q CD-like SU(iV c ) gauge theories with Nf flavors of quarks Q and antiquarks Q, any 
charges (like baryon number) associated with conserved vector currents are not sponta- 
neously broken Thus, the lightest bound states of Q's or Q's carrying nonvanishing 
baryon numbers are almost stable and they will decay into the ordinary quarks and lep- 
tons through some baryon-number violating nonrenormalizable operators suppressed 
by the gravitational scale M* ~ 2.4 x 10 18 GeV. If the dynamical scale A of the SU(iV c ) 
gauge interactions is well below the gravitational scale and N c is sufficiently large, the 
lifetimes of the superheavy baryons may be longer than the age of the present uni- 
verse. Therefore, the presence of long-lived superheavy baryons is a generic prediction 
of QCD-like gauge theories for a certain parameter region of iV c and A. 

However, any (quasi-) stable particles much heavier than 0(1) TeV are cosmo log- 
ically dangerous, since they would easily overdose the universe if they were once in 
the thermal equilibrium ||. One may usually invoke some inflationary stage in the 
universe's evolution to dilute the number density of such superheavy X particles. If 
the reheating temperature after the inflation is much lower than the masses of X par- 
ticles one may completely neglect the thermal production of X particles. It has been, 
however, suggested || [| that if the X-particle masses mx are of order of the Hubble 
constant H at the final epoch of inflation, gravitational interactions may give nonneg- 
ligible contributions to the X-particle production just after the end of the inflation. 
The numerical calculation in Ref. 0, in fact, shows that when mx — (0.04 — 2) x H 
a suitable amount of the X particles is produced to form a part of the dark matter in 
the present universe. It is remarkable that the decays of such X particles will generate 
significant effects on the spectrum of high energy cosmic ray if the lifetimes of the X 
particles are of order of the age of the present universe. 

Several events of the ultra-high energy (UHE) cosmic rays || ^, [7] have been recently 
observed beyond the Greisen-Zatsepin-Kuzmin (GZK) bound E ~ 5 x 10 10 GeV ||. 
These are naturally explained ||, |1(| by the decay products of superheavy X particle 
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of mass m x — 10 13 — 10 14 GeVQ with lifetime r x — 10 10 — 10 22 years] 2 ] if its energy 
density px lies in the range 

Px / Pc ~ io- 12 - 1. (1) 

Here, p c ~ 8.lh 2 x 10~ 47 GeV 4 with h ~ 0.5 — 1.0 is the critical density of the present 
universe. Since the required window of the energy density px is very wide, this scenario 
seems very plausible and attractive. 

In this paper, we construct explicit models based on supersymmetric (SUSY) SU(iV c ) 
gauge theories with Nf flavors of quarks Q and antiquarks Q, in which the lightest 
baryons B and antibaryons B have the desired mass and lifetime, that is, mg = ~ 
10 13 — 10 14 GeV and tb = t b ~ 10 10 — 10 22 years. In these models the gaugino con- 
densation necessarily occurs and hence the models may play a role of so-called hidden 
sector for SUSY breaking in supergravity [13]. Thus, the long-lived superheavy B and 



B are regarded as byproducts^ of the hidden sector gauge theories for dynamical SUSY 
breaking.^] 

It should be noted here that in contrast to the non-SUSY case, the SUSY QCD-like 
gauge theories may yield baryon-number violating vacua due to the presence of scalar 
quarks. In these vacua we have no longer quasi-stable baryons. However, if quarks Q 
and antiquarks Q have SUSY-invariant masses, the unwanted baryon-number violating 



vacua disappear as shown in Ref. [16|. Thus, we consider the SUSY QCD-like gauge 
theries with massive quarks throughout this paper. 

In section 2, we briefly discuss vacua of SUSY SU(iV c ) gauge theories with Nf pairs 
of massive quarks and antiquarks. We restrict our discussion to the case of Nf = N c + 1 
and show that there is a unique SUSY-invariant vacuum preserving the baryon-number 
conservation. Thus, we always have stable baryons and antibaryons in this theory. In 



1 There is an analysis which suggests mx — 10 12 GeV [O]. 

2 The required lifetime may be accounted for by imposing discrete gauge symmetries even if the 
X are elementary particles jL2j . 

3 A similar idea has been considered in connection with string theory |l4| . 

4 We may consider non-SUSY SU(iV c ) gauge theories which cause dynamical breaking of the Peccei- 
Quinn symmetry at A ~ 10 13 GeV Jig ]. In these models we may have naturally quasi-stable baryons 
of masses ~ 10 13 — 10 14 GeV, which have the required lifetimes tb — 10 10 — 10 22 years for N c ~ 5, 6. 
The main decay mode of such baryons will be B — »• Z+Higgs boson or 2xHiggs bosons. 
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section 3, we extend the above model to the supergravity, in which we introduce non- 
renormalizable operators. We show that possible baryon-number violating operators 
induce spontaneous breakdown of the baryon-number conservation and the baryon- 
meson mixings occur. Owing to the baryon-number violating effects even the lightest 
baryons are no longer stable. However, we see that the lifetimes of the baryons can be 
chosen as in the required range to account for the UHE cosmic rays by taking N c = 8, 9 
and 10. We also show that the gaugino condensation necessarily occurs and the SUSY 
may be broken in the dilaton stabilized vacua. In section 4, we argue that the lifetimes 
of the baryons become longer if the baryons carry nonvanishing charges of some extra 
symmetries. As for such symmetries we adopt the matter parity Z2 or the discrete 
baryon parity Z 3 , since they are often used to guarantee the stability of usual proton 
|T7| 1. In these cases we find the desired lifetimes are obtained for somewhat smaller 



N c = 6 — 9. The last section 5 is devoted to discussion and conclusions. 

2 Supersymmetric QCD with Massive Quarks 

Let us consider SUSY SU(iV c ) gauge theories with Nf flavors of quarks Q % a and an- 
tiquarks Q%, where a = 1,---,N C and i,i — 1, * • • , Nf. We omit the color index a, 
hereafter. We neglect the mass term for Q l and Qi for the time being. Then, we have 
a global SU(iV/)xSU(iV/)xU(l)yxU(l)ij symmetry. We restrict our consideration to 
the case Nf — N c + 1, since the dynamics is the clearest in this case. 

For Nf = N c + 1, the low energy physics is described by canonically- normalized 
gauge invariant composite fields, mesons Ml ~ Q t Qj/A, baryons Bi ~ eijk-.-iQ^Q k • • • Q'/A^ -1 
and antibaryons B l ~ e^ k '" l QjQ% ■ ■ ■Qi/A Nc ~ 1 ||18|| . The dynamically generated super- 
potential is given by 

W dyn = BiM}B' - ^^detM. (2) 

The mesons M|, baryons Bi and antibaryons B l are all massless to satisfy the 't Hooft 
anomaly matching conditions [18[ . 



Now we introduce the mass term for quarks Q l and antiquarks Q^. Then, the total 
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effective superpotential is given by 



W = B t M l j& - ^^detM + m\kM). (3) 



It is a straightforward task to see a SUSY invariant vacuum: 



(Ml) = AT(detm)^(m- 1 );-, (4) 
(Bt) = iff) = 0. (5) 



In this vacuum the mesons Mj, baryons Bi and antibaryons B l have the following mass 
terms 



W mass = A^detm^m- 1 )^^" 



l-A^c (detm) nI {m\m k - m 3 k m l ^)M\Mi . (6) 
2 J 



Thus, we obtain the masses for these composite fields as 



m B ~ (mh Nc ~ 1 )^, (7) 
m M ~ (m C-1 A) : 



i 



where we have assumed a common mass m:- = m5f, for simplicity. We will identify 
these composite baryons Bi and antibaryons B l with the long-lived superheavy X 
particle introduced to explain the UHE cosmic rays discovered beyond the GZK bound. 
Therefore, we take 

m B ~ (mA^- 1 )^ ~ 10 13 - 10 14 GeV. (9) 
Using the Konishi anomaly relation jl9] we determine the gaugino condensation as 



(AA) = ^ml^Qj) = A^(detm)^ 

(10) 

~ (m iVc+1 A 27Vc ~ 1 )^. 

This condensation may give a dominant contribution to the SUSY breaking in super- 
gravity. 
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3 Extension to the Supergravity 



We now extend the above model to the supergravity. So far we have considered only 
renormalizable interactions, but in the framework of supergravity it is quite natural to 
introduce nonrenormalizable interactions. Namely, we introduce in general^ 

W tIee = mjQ'Qj + 1 ^- E e ijk ... l Q j Q k ■■■Q l + -j-^e^QjQ, ■ ■ ■ Q h (11) 

with bi,bi = 0(1). We have omitted the other nonrenormalizable terms which are 
irrelevant for our purposes here. Then, the total effective superpotential is given by 

w = w dyn + w tiee 

+m ; AM j + 6i (_) A^, + (*(-) KB 1 . (12) 



It is a straightforward task to see a SUSY invariant vacuum: 



(Mj) = A^(detm)^(m- 1 );-, (13) 
\M, 



(Bi) = -iA^-) A*(detm)^mf, (14) 



— \ /A \ Nc ~ 3 v c+ i -i _ 

B l ) = -V j A— (detm)^m}, (15) 

up to the leading order in A/M*. In this vacuum the mesons Mj, baryons Bi and 
antibaryons B % have the following mass terms 

-~A^(detm)^(m^ml - m\m l ^)M\Mj 
2 J 



, A siVc-3 

-b k { — ) A— (detm^miBiM'j 



KM, 

( A \^- 3 m 



AT. *n 



-b- k [ — ) A— (detm)^mfM^". (16) 



There exist also baryon-number violating nonrenormalizable operators in the Kahler poten- 
tial. However, they are negligible compared with the baryon-number violating operators in the 
superpotential. 
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Thus, we obtain mixing masses between the composite meson and baryon fields as 

A \ Nc ~ 3 



m BM ~ ^-j (m- 1 A^ +1 )^. (17) 

The diagonal masses tub for baryons and tum for mesons are the same as in Eq. (§) up 
to the leading order in A/M*. Notice that as long as A -C M* and N c ^> 3 the mixings 
between mesons and baryons are very small.Q 

Let us now discuss the decay of these Bi and B l . When the masses ttim for the 
mesons are all larger than the half of those of Bi and B l , these composite baryons should 
decay directly into ordinary light particles including quarks and leptons through the 
mixing terms in Eq. ([TED together with the following nonrenormalizable operator: 

W = -^QQHH. (18) 

The lifetimes of Bi and B % are determined as 

1 / M \ 2 (^- 2 ) ! 

1.4t (m- 3 A^ +3 )^, (19) 



/2 V A 

which should be taken ~ 10 10 — 10 22 years to account for the UHE cosmic rays. 

When 2wm < tub, we have new decay channels Bi (B l ) — > 2M|. Since the mesons 
Mj decay into ordinary light particles very quickly through the interactions Eq. (|18"D, 
the lifetimes of Bi and B % decaying into the ordinary light particles are determined by 
the 2Mj decay channels which are given by 

, M x2(V c -3) 

* N ' ,™ c -7 \-N c +7\j£ 



TB ~ ( ^ ) (m 2 ^- 7 A-^ +7 ) . (20) 



These are somewhat shorter than the previous lifetimes Eq. (|19f). Thus, we conclude 



' B 



10-54 _ 10 -42 GeV ( for 2m M > mB ) 

( ] A_) 2(7Vc " 3) (m 2iY c -7 A -V c+7) ^ 

1(T 54 - 1(T 42 GeV (for 2m M < m B ) 



(21) 



(> 



The baryon-number condensations in Eqs. ( J14| , |15| ) induce kinetic mixings between the meson 
and baryon fields in the Kahler potential, which give the same-order effects as those discussed in the 
text. We neglect them, for simplicity, since they do not affect our main conclusions. 
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to have the required lifetime r# ~ 10 10 — 10 22 years. 

Now we are at the point to discuss the gaugino condensation Eq. ([10|) . It is well 
known that this gaugino condensation may induce SUSY breaking together with the 
dilaton field stabilization in supergravity (T^J.Q Assuming the gravitino mass m 3 / 2 — 



1 TeV we determine the gaugino condensation scale (AA) 1 / 3 ~ 10 13 GeV. From the 
constraints Eqs. (|S|, ID, |2T[) we obtain the following relations: 



2iv2-iY c +6 



^M^-^AA)^^ 1 = m B ^" 2 (fov2m M >m B ), 

9 ( 22 I 

2 2JV^-9iV e + 16 V I 

Ml {Nc ~ z) {\\)—t b 1 = m B Nc ~ 2 (for2m M <m B ). 

These relations are consistent with the required values of rriB, t b and (AA) only when 
the numbers of colors N c are 

N c = 8,9,10 (for 2m M > m B ), 

N c = 9,10 (for2m M <m B ). 1 ' 

Here, we have assumed the coupling constant f — 1. Then, A and m are determined 

as 

A ~ 10 13 - - 10 14 - 5 GeV, (24) 
m ~ 10 10 - 5 - 10 13 - GeV. (25) 

Note that these numerical values yield too large a mass term for the Higgs doublets in 
Eq. (|T8"D. We postpone the discussion on this point to the final section. 

4 Models with Discrete Gauge Symmetries 

In the previous section we find that the desired lifetimes of the baryons are obtained 
if N c = 8, 9 and 10 without invoking any extra symmetries. In this section, we impose 
the matter parity Z2 or the baryon parity Z3 on our model, since these discrete gauge 
symmetries are often used to suppress very rapid decays of the usual proton.^ The 
charges for the minimal SUSY standard-model (MSSM) particles under the matter 
parity Z2 and the baryon parity Z3 are given in Table 0. If the composite baryons do 

7 The SUSY-breaking effects do not change the order of magnitude of (AA) in Eq. (|l~0|). 

8 In the case that neutrinos acquire Majorana masses through operators W = {1/Mr)UHH ', the 
anomaly-free discrete gauge symmetries are only the matter parity Z2 and the baryon parity Z3 ]12fl . 
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z 2 
z 3 



q u d I e H H 

o I I o i i r 

2 1 2 2 1 2 



Table 1: Charges for the MSSM particles under the matter parity Z 2 and the baryon 
parity Z 3 . q, u, d, I and e denote SU (2) L-doub\et quark, up-type antiquark, down-type 
antiquark, SU(2) ^-doublet lepton and charged antilepton chiral multiplets. H and H 
are chiral multiplets for Higgs doublets. 

not carry nonvanishing charges of the Z 2 and Z 3 , the analyses are the same as in the 
previous section. If these composite baryons have nontrivial Z 2 or Z3 charges, however, 
there are no mixing mass terms between mesons and baryons, since the linear terms 
of Bi and B % in Eq. ( ff2|) are forbidden {i.e. b % = b l = 0) and hence (Bj) = (S 1 ) = 0. 
Thus, the decay channels of Bi and B l are different from those in the previous section. 



In this section, we discuss the lifetimes of the baryons with the discrete Z 2 or Z3 fllzfl . 

First, we consider the case of the matter parity Z 2 . The charges for the Bi and 
B l must be opposite for quarks Q l and antiquarks Qj to have invariant masses. Thus, 
we suppose that both Bi and B l are odd under the Z 2 . Then, the lowest dimensional 
operators which cause decays of the composite baryons are 

W = j^QQ ■■■QIH+ j^QQ ■■■QlH, (26) 

where Q's and Q's, for example, carry Z 2 charges (2r + l)/iV c (r G Z) and — (2r + l)/N c , 
respectively.^ Then, the lifetimes of Bi and B % are determined as 

tb * (jf) H^ 1 )^. (27) 

Next, let us turn to the case of the baryon parity Z3. We suppose that the Bi carry 
+ 1 and the B l carry —1 of the Z 3 charges. Then, the lightest baryons decay into the 
MSSM particles through the following operator: 

W = j^ K QQ---Qudd, (28) 



The Z2 is anomaly free. 
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where Q's and Q's carry Z 3 charges (3r + 1)/N C and — (3r + 1)/N C , respectively.^ The 
lifetimes of 5j and 5* are given by 

/ A/7" \ 2-/V c 

r B ~(^J (m 3 A^" 3 )^. (29) 

Note that even if we assign the Z 3 charge —1 for B^ and +1 for B l , the lightest baryons 
decay through the operator 

W = j^ K QQ---Qudd, (30) 

so that the lifetimes are the same as in Eq. (|2"P|). 

From the constraints Eqs. (§, 0, |27|, |29|) we obtain the following relations: 

2Jv2+JV c -4 

(f^y 7. A 

(31) 



o^)\r n 2(JV c -l) . 

M^-V (M)-^ = m B N ^ 2 (forZ 2 ), 



2(JV C -1) 



M ™° (XX)^^ t b 1 = m B (for Z 3 ), 

which are consisitent with the required values of m^, tb and (AA) only when the 
numbers of colors N c are 



N c = 7,8,9, (forZ 2 ), 

iV c = 6,7,8 (forZ 3 ). 1 ' 

Therefore, when the Bi and B l carry nonvanishing charges of the matter parity Z 2 or 
the baryon parity Z3, the required numbers of colors iV c to obtain the desired lifetimes 
are smaller than those in the case without Z 2 or Z 3 . Notice that A and m are almost 
the same as in Eqs. (|4], ^5)) inspite of the change of iV c . 



5 Discussion and Conclusions 

In this paper we have constructed explicit models based on SUSY SU(iV c ) gauge the- 
ories with Nf flavors of quarks and antiquarks which naturally accommodate the su- 
perheavy composite baryons B and B introduced to account for the UHE cosmic rays 
beyond the GZK bound. The models contain three crucial parameters, the quark mass 
m, the dynamical scale A and the number of colors N c . The number of flavors Nf 
is fixed as Nf = N c + 1, for simplicity. In these models the gaugino condensation 
10 The Z3 is anomaly free. 
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(AA) always occurs, which may cause the SUSY breaking in supergravity. Assuming 
(AA) 1/3 ~ 10 13 GeV (for m 3 / 2 — 1 TeV) and the desired properties for the composite 
baryons (mg = m§ ~ 10 13 — 10 14 GeV and tb = t b ~ 10 10 — 10 22 years), we have 
obtained m ~ 10 10 5 - 10 130 GeV, A ~ 10 13 - 10 14 - 5 GeV and N c = 6 - 10. Namely, 
we have found that the required long lifetimes of superheavy B and B are naturally 
explained in SUSY QCD-like gauge theories with large number N c of color degrees of 
freedom. Although we have restricted our analysis only to the case of Nj — N c + 1, it 
is possible to consider other cases. 

Several comments are in order. We should mention first the so-called [i term prob- 
lem. Owing to the QQ condensation, the Higgs doublets seem to have an invariant 
mass term fiHH, where /J ~ / (QQ^ /M# ~ (10 8 — 10 10 ) x / GeV. Thus, in order to 
induce a correct vacuum-expectation values for the Higgs doublets, a negative invari- 
ant mass of order 10 8 — 10 10 GeV must be introduced to cancell the unwanted large 
mass / (qQ^ /M*. Alternatively, one can solve this problem by putting the coupling 
constant / of the operator (1/M*)QQHH very small / ~ 10~ 5 — 10~ 7 . We see that 
even if it is the case, the obtained A and m are almost the same as in Eqs. (|24| , p5|) , 
and our conclusion does not change much.[^] The small value of / will be understood 
by some axial symmetries, which may also explain the small mass m for Q l and 

Our model is also applicable to the gauge- mediated SUSY breaking scenario |2(J . If 
the gaugino condensation causes SUSY breaking partially, the induced gravitino mass 
"7713/2" must be smaller than 1 GeV in order to suppress dangerous flavor-changing 
neutral currents sufficiently We find that the desired composite baryons are obtained 
if m ~ 10 2 - - 10 9 - 7 GeV, A ~ 10 13 ' 3 - 10 16 - GeV and A^ c = 6 - 11 for "m 3/2 "~ 

11 For the case of 2tom > ttib, if we take / ~ 10~ 5 — 10~ 7 the dominant operators contributing to 
the baryon decays are not those in Eq. ( fj"8| ) but direct decay operators given by 

W=—^- T QQ---QHH+—j^ T QQ---QHH, (33) 

which give the desired number of colors N c = 7, 8, 9 instead of 7V C = 8, 9, 10 obtained in the text. On 
the other hand, the lifetimes are independent of the value of / for the case of 2tom < tub- 

12 R symmetry may be an example in which HH has i?-charge zero. We naturally obtain the 
terms mQQ(l + f'HH/M^) in the superpotential with a coupling /' of order one, which may yield an 
appropriate fj, term. The small mass m is regarded as a breaking term of the R symmetry. To have 
unsuppressed baryon-number violating operators in Eq. (|33|), for example, we assume that both of Q l 
and Qi carry i?-charge 2/N c . 
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100 keV - 1 GeV. 

Finally, we should comment on the gauge coupling constant a c of the SUSY SU(iV c ) 
gauge theories considered in this paper. With N c = 6 — 10 and Nj = N c +1 the solution 
to the one-loop renormalization group equation for the gauge coupling constant a c is 
given by 

a -i(A) =a -i(M.) + ^ln(A). (34) 

Using a c (A) ~ oo we get a c (M*) ~ 0.03 — 0.05 at the gravitational scale M*. It 
is interesting that the value of a c (M*) is roughly consistent with the hypothesis of 
unification with the standard-model gauge coupling constants. That is, the larger 
dynamical scale A ~ 10 13 — 10 14 GeV compared with the usual QCD scale A QCD ~ 
0.1 GeV is a natural consequence of the large number of colors N c = 6 — 10. 
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